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Abstract 

We seek exact solutions to the Einstein field equations which arise when two spacetime 
geometries are conformally related. Whilst this is a simple method to generate new 
solutions to the field equations, very few such examples have been found in practice. We 
use the method of Lie analysis of differential equations to obtain new group invariant 
solutions to conformally related Petrov type D spacetimes. Four cases arise depending on 
the nature of the Lie symmetry generator. In three cases we are in a position to solve the 
master field equation in terms of elementary functions. In the fourth case special solutions 
in terms of Bessel functions are obtained. These solutions contain known models as special 
cases. 



1 Introduction 



Exact solutions, as opposed to computer generated solutions, of the Einstein field equations are 
of immense importance in understanding the behaviour of a large variety of celestial phenomena. 
The literature abounds with many different techniques that have been invoked in an effort to 
obtain new exact solutions for different configurations of matter. 

Conformal transformations, as a mathematical procedure, have also been successfully utilised 
in obtaining new solutions. This is amply illustrated by the comprehensive model of Castejon- 
Amenedo and Coley (1992). Further, it is known that conformal structures play an important 
role in twistor theory (Penrose 1999). The restrictive feature of such analyses, however, is the 
complexity and nonlinearity of the resultant field equations. Some researchers have adopted 
the Newman-Penrose formalism which, on account of severe integrability constraints arising 
out of the Bianchi identities, have not proved fruitful in general. 

An alternative is to assume the existence of a conformal symmetry on the manifold in order 
to generate solutions to the Einstein field equations, as these symmetries impose additional 
restrictions on the metric tensor field and the field equations may be simplified. The physical 
significance of conformal Killing vectors is that they generate constants of the motion along 
null geodesies for massless particles. 

The maximal spanning G15 of conformal motions for Minkowski space is given by Choquet- 
Bruhat et al (1982), and for Robertson- Walker spacetimes by Maartens and Maharaj (1986). 
In addition the conformal geometries of the pjo-wave spacetimes (Maartens and Maharaj 1991), 
static spherically symmetric spacetimes (Maartens et al 1995,1996), Bianchi I and V locally 
rotationally symmetric spacetimes (Moodley 1992) and Stephani spacetimes (Moopanar 1993) 
have been completely determined. Kramer (1990) was able to generate a class of metrics for 
rigidly rotating perfect fluids which admit a proper conformal Killing vector. Additionally 
rigidly rotating perfect fluids admitting two Killing vectors and a proper conformal Killing 
vector were studied by Kramer and Carot (1991). Axisymmetric spacetimes, in the general 
case of differential rotation, were examined by Mars and Senovilla (1993, 1994). 

Several new results, utilising this conformal symmetry approach have recently been obtained 
by Castejon-Amenedo and Coley (1992), Coley and Tupper (1990a,b,c), Dyer et al (1987) 
and Maharaj et al (1991). In particular, Maartens and Mellin (1996) used the conformal 
symmetries of Bianchi I spacetime to demonstrate that the expansion of anisotropic radiation 
universes tends towards isotropy at late times. 

Our interest lies in spacetimes that admit an s-dimensional Lie algebra C s of conformal motions. 
We exploit the Defrise-Carter theorem (1975) to generate new models of perfect fluids. We 
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select a spacetime of Petrov type D and its conformally related counterpart in an effort to 
obtain new solutions to the associated Einstein field equations. We analyse the field equations 
and show that they can be reduced to a simpler form. The results of Castejon-Amenedo and 
Coley are regained as a special case of a more general class of exact solutions. Group invariant 
solutions to a particular field equation which acts as a master equation for the entire system are 
sought. This analysis generates a rich class of solutions. A number of cases arise depending on 
the nature of the Lie symmetry generator. In all cases we are in a position to provide solutions 
to the master field equation in terms of elementary functions. This analysis demonstrates the 
value of the Lie analysis of differential equations in this application. 



2 Spacetime Geometry 

We consider the line element 

ds 2 = -dt 2 + dx 2 + e 2uM (dy 2 + dz 2 ) (1) 
which is of Petrov type D. The spacetime (JIJ admits three Killing vectors 

d d 

Xs = x Ft +t 3- x (4) 

which obey the commutation relations 

[Xi,X 2 ] = [X!,X 3 ]=0 [X 2 ,X 3 ]=0 

The Lie algebra of Killing vectors is a G3 of motions, with the group structure satisfying the 
above relations. 

Suppose that a manifold (M, g) is neither conformally flat nor conformally related to a gener- 
alised plane wave. Then according to a theorem due to Defrise-Carter (1975), a Lie algebra of 
conformal Killing vectors on M with respect to g can be regarded as a Lie algebra of Killing 
vectors with regard to some metric on M conformally related to g. Therefore if a spacetime ad- 
mits the conformal group C s , then either it is conformally flat (s = 15), conformally related to 
a generalised plane wave (s < 7), or the metric g a b = e 2U g a b where g a b admits an s-dimensional 
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(s < 6) isometry group. The last possibility is evident above and so as our starting point we 
consider the conformally related metric 



ds 2 = e 2U 



. d 1? + dx 2 + e 2v ^ (dy 2 + dz 2 ) 



(5) 



where U = U(t,x,y,z). The Killing vectors (J2J) are now conformal Killing vectors of the 
conformally related spacetime ©. The Weyl conformal tensor is given by 



3e 2t/ C, 



0101 



-6C, 



0202 



-6C, 



0303 



6C 



1212 



— 6C*i3i3 — —e 2l/ C2323 — -6 2i/ C3232 — ^j/j/ + ^zz 

which clearly indicates that the metric (jSJ) is not conformally flat in general. (We use the 
notation that overhead bars on quantities are defined in the conformally related spacetimes 
(jHJ).) Note that v m + v zz is not zero. We make the assumption that 

Vyy + Vzz = -2ke 2u (6) 

where k is a nonzero constant. This choice is made on the grounds of simplicity and follows the 
treatment of Castejon-Amenedo and Coley (1992). In addition if k = then the line element 
(JTJ) becomes flat, and consequently (J3J would be conformally flat. The condition (jHJ), with 
k 7^ 0, obviates this occurrence. 

To determine the perfect fluid energy-momentum tensor, we select a fluid 4-velocity vector u 
that is noncomoving with the form 

u a = e~ u (cosh v5% + sinh vSf) (7) 

where v = v(t, x). Note that a trivial calculation reveals that for v = constant, the perfect fluid 
Einstein field equations would imply conformal flatness. 



3 Field Equations and Kinematics 

The Einstein field equations are given by 

U t U y -U ty = 



U t U z -U tz = (9) 



3 



U X U Z - U s 







u t u x - u tx 



(fj, + p)e sinh2t> 



UyU z 



Uy Z +U z Uy + VyU Z 







-2U XX - U 2 + W 2 



-e 2v (2U yy + 2U ZZ + U 2 + U 2 + Vyy + v = (/j, + p)e 2U cosh 2 v - pe 



+e' 2u (2U yy + 2U ZZ + U 2 + U 2 + Vyy + v ' M ) = (fi + p) e 2U sinh 2 v + pe 



2U ZZ + U 2 + 3U 2 + 2u y U y - 
2v z U z + e 2v (2U XX - 2U tt + U 2 X - U 2 ) = pe 

2Uyy + Ul + 3U 2 ~ 2VyUy + 2V Z U Z 

+e 2v (2U XX - 2U tt + U 2 - U 2 ) = pe 



so that the gravitational field is irrotational. The components of the acceleration ii 



2U tt - U 2 + 3U 2 



for the line element 



The vorticity uj a b is given by 




Ul = 



(v t + U x ) cosh 2 v + (v x + U t ) cosh v sinh v 



which are nonzero in general. The expansion is given by the expression 

6 = e~ u [(v t + U x ) sinhw + (v x + U t ) coshf] 

The shear tensor components have the form 
1 



0"n 



(2 sinh 2 v - l)e u + e~ u [(v t + U x ) sinh v + (v x + U t ) cosh v] 



a i = — e u sinh v cosh v [(v t + U x ) sinh v + (v x + U t ) cosh v] 
3 



i 



(2 cosh 2 v + l)e u - e~ u [(v t + U x ) sinh v + (v x + U t ) cosh v] 



(20) 

(21) 
(22) 
(23) 



0-22 = ~^f V U K v t + U x) sinh v + (v x + U t ) cosh v] 



^33 



(24) 



which do not vanish in general. From p8|) - (|2*T|) we observe that the fluid congruences of the 
conformally related line element (0) are accelerating, expanding and shearing. 



4 Reduction of the Field Equations 



The most general functional form admitted by the equations is 

e~ u = f(t,x) + h(y,z) 



(25) 



where / and h are arbitrary. The remaining field equations (|12)l - (|17|l . respectively, then assume 
the following form 

2(f + h)ftx = — (/i + p) cosh f sinh v (26) 



<L)Z 



V z ky + Uyk Z 



(27) 



"3 (f? ~ fl) ~ 2(/ + h)f„ - 2k(f + hf 



+ e 



-2v 



3(h 2 y + hl)-2(f + h)(h yy + h zz )] = -(fx + p)cosh 2 v+p (28) 



-3 (/ 4 2 - fl) + 2(/ + h)f u - 2k(f + hf 
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+ e 



-2v 



3 (hi + hfj - 2{f + h) {hyy + h ZZ ) 



[jj, + p) sinh 2 i> + p 



(29) 



+ e 



-3 (f? - /*) + 2(/ + (/« - /„) 
3 (/ij + ^) - 2(f + /i) - v z h z + 

"3 (ft - fl) + 2(/ + h) (ftt - f„) 



P 



+ e 



-2v 



3 (^ + K) - 2(f + h) (v z h z - v y h y + h 



y u y 1 To/ 



P 



where the variable U has been replaced by / and h via (|25|) . 



(30) 



(31) 



The dynamical quantities have the following forms: 



P = -3(f?-f*)+2(f + h)(f tt -f x 



+ e~ 2v (3 (/i 2 + hi) ~(f + h) (hyy + h z , 
for the pressure p and 

A* = 3 (/ 2 - / 2 ) + 4fc(/ + /i) 2 - 3e- 2 ^ + hi) - (/ + &)(/»„, + ft, 



(32) 



(33) 



for the energy density /i. We now seek an expression for the quantity v(t,x). This is accom- 
plished by first substituting in (J28)l for p (but we retain the p cosh 2 t> term) to give 



(ji + p) cosh^ v = 2(f + h)f u + 2k(f + h) 2 + eT lv (f + h) (h yy + /i ; 
On substituting (J32j) into (|29j) for p (but we retain the p sinh 2 v term) we get 
(ji + p) sinh 2 v = 2(f + h)f xx - 2k(f + h) 2 - e- 2v (f + h) (h yy + h 



Dividing fl35J) by (jHU) yields 



tanh 2 v 



2f xx - 2k(f + h)- e- 2v (h V y + h 



2 ftt + 2k(f + h) + e- 2u (h yy + h zz ) 
Another expression for tanh v may be obtained by dividing (J33j) with (j2T?j) to give 

(2f xx - 2k(f + h)- e- 2v (hyy + h zz )f 



tanh v 



4/ t 2 



Comparing this equation with (|3l)j) yields the differential equation 

f 2 x = (2f xx -2k(f + h)-e- 2 »(h yy + h zz : 
x (2/« + 2fc(/ + /i) + e- 2 "(/i ra + h z 



(34) 



(35) 



(36) 



(37) 



On rearranging (|36|) we obtain 

2fxx — 2ftt tanh 2 v -iv 11, , u \ 

2k t = 2kh + e h vv + n z , 

tanh 2 z; + l ^ v vy ' 

from which it is clear that the left hand side is a function of t and x whereas the right hand 
side is a function of y and z. This implies that the variables separate and we can put 



2/ CT -2/ tt taah 2 i> = ^ + & _ 2v 

tanh u + 1 

= a 

where a is a constant. We thus derive the expression 

h yy + h zz = e 2u {a - 2kh) (38) 
Equation (|HK|) further simplifies the system pK |) -(|3*T ]) . 

From the above analysis it is clear that the field equations (j2SJ)-(JSIJ) can be expressed in a 
simpler form. The resulting system is given by 

V = 3(f?-f 2 x ) + (f + h)(4kf-2kh + 3a)-3e- 2 »(h 2 y + hl) (39) 

p = -3(f 2 ~f 2 ) + (f + h)(2f tt -2f xx + 2kh-a) 

+ 3e- 2 " (h 2 y + h 2 z ) (40) 

tanh 2 , = 2/ f ~ 2 ; fc { ~ a (41) 
2f tt + 2kf + a 1 ' 

ftx = \(2f xx -2kf-a)(2f tt + 2kf + a) (42) 
h yz = v z h y + v y h z (43) 

hyy - h zz = 2v y hy + 2v z h z (44) 

The system subject to condition (jUJ), viz. v yy + v zz = —2ke 2u , must be solved in 

order to generate a conformally related perfect fluid model. 

Castejon-Amenedo and Coley have considered the case h = constant; this constant may be 
effectively absorbed into / without any loss of generality and we can consequently set 

h = 
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With this value of h the Einstein field equations reduce to 

P = -3 {ft ~ fl) + Viftt ~ /„) (45) 

V = Z{fi-fl)+±kf 2 (46) 

v = t a nh- 1 ( kf ~ fxx ) (47) 

\ Jtx J 

fl = (fu + kf) (/„ - kf) (48) 

which is in agreement with the equations used by Castejon-Amenedo and Coley (1992). Note 

that corresponds to h — a — 0. It is possible in P^ - pijl to have a = with h ^ 0. 
This naturally leads to two categories of exact solutions that we now present. 



5 An Extension of the Castejon-Amenedo and Coley so- 
lutions: h = 

This category of solutions corresponds to h — a = 0. In order to generate a solution to the 
system fl43 j) - PHj) it is sufficient to obtain a form for / = f(t,x); in other words (|4*Hj) must be 
integrated. To obtain a solution Castejon-Amenedo and Coley have assumed that the variables 
separate and set 

f(t,x) = F(t)G(x) (49) 

which yields the differential equations 

FF — ■yF 2 + kF 2 = (50) 

GG" - -G' 2 - kG 2 = 0, (51) 

7 

where 7 is a constant. 
5.1 CASE I: 7 = 1 

It is easy to solve equations (j5Uj) - (j5Tj) in general. As a result, the system (j43|) - (j4Tj) reduces to 

~3k 2 (t 2 - x 2 ) - 6fc(fcit + k 2 x) + 3(k 2 - k 2 2 ) + Ak 



^ = e k(x 2 -t 2 )+2k 1 t+2k 2 x 



(52) 
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p 



^k(x 2 -t 2 )+2k 1 t+2k 2 x 



-k 2 {t 2 - x 2 ) + 2k(ht + k 2 x) - (kf - k 2 ) - 4k 



(53) 



tanh v 



kx + k 2 



kt — k\ 

This general solution contains the solution of Castejon-Amenedo and Coley when we set 

k x = k 2 = 

in (|52j) - (|54j) . For this choice of constants we obtain 



(54) 



D k(x 2 ~t 2 ) 



3k 2 (t 2 -x 2 ) +4k 



(55) 



D k{x 2 -t 2 ) 



k 2 (x 2 - t 2 ) - 4k 



(56) 



x 

tanhf = — (57) 

which is the exact solution of Castejon-Amenedo and Coley (1992). Consequently we have 
extended their solution (|55 j) -(|57 }) to the more general class (|52)l - (j54~jl . 

For the exact solution (|52j) - (|54j) we obtain 

fi + p = 2e k ^ 2 - t2 ^ t+2k ^\k\t 2 -x 2 )-2k(k 1 t + hx) + (kl-kl 



_Q] ie k(x 2 -t 2 )+2k 1 t+2k 2 x 



[i-p = 4 :e k(x 2 ~t 2 )+2k 1 t+2k 2 x 



~k 2 {t 2 - x 2 ) - 2k{kit + k 2 x) + (k{ - k 2 2 ) + 2k 

and it is possible to study the weak, dominant and strong energy conditions. We note that 
the appearance of the constants kx and k 2 allows for a wider range of behaviour for our class 
of solutions than is the case with the Castejon-Amenedo and Coley (1992) exact solution. For 
this model u a = so that the field is nonaccelerating; however both a a t and G are nonzero 
which implies that the gravitational field is shearing and expanding. 



5.2 CASE II: 7 ^ 1 

The case with 7 7^ 1 in (J45|) is of interest in generating cosmological models for investigating a 
wider range of physical behaviour. The transformations 

F = wp, G = uf T 
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reduce equations (J50*j) - (J5T|) to 

iii + k(l — j)ui = 



tin + fc— — —u 2 = 

7 

which are linear in u\ and 1x2 respectively. The solutions of this system of differential equations 
are easily obtained in terms of elementary functions and consequently we do not present their 
explicit forms here. 



6 Group Invariant Solutions: h ^ 

This category of solutions corresponds to h ^ with a = 0. It is possible to solve the more 
general system (|39j) - (|4"4"j) with h ^ 0. If we take h = h(y), then (|43|) implies that v = v{y)- 
Consequently, the condition (JHJ) becomes the ordinary differential equation 

v yy = -2ke 2u . (58) 

Equation (}4"4"j) reduces to 

h y = Ce 2u , (59) 
where C is a constant of integration. The general solution to (|HHj) and f!59j) is given by 

_ __ 

A; ' 



( - = - -cosech 2 (y/2Ay + b) (60) 



Ky) = ^^coth(V2Ay + B)+D, (61) 

where A, B and are constants of integration. In order to complete the solution we need to 
solve (J42j) and obtain / as this immediately yields expressions for p, fj, and v. We now focus 
our attention on (|4"2jl with a = 0, i.e. 

fi = (fu + kf) (f xx - kf) (62) 

with the objective of generating solutions in a systematic manner. We achieve this by invoking 
the method of Lie group analysis of differential equations (Olver 1993). The Lie symmetry 
generators of (J62)) are found with the help of the computer package PROGRAM LIE (Head 1993). 
The infinitesimal Lie symmetry generators for the partial differential equation ()62)1 are 

Z. = I (63) 
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^ = f§j (66) 

It is clear that these symmetries form the Lie algebra is A^@A\ (Patera and Winternitz 1977) 
with basis given by 

Gi = Zx + Z 2 (67) 

G 2 = Z x -Z 2 (68) 

G 3 = Z 3 (69) 

G 4 = Z 4 . (70) 

It is this basis that we will use in our subsequent analysis. Before proceeding further, we observe 
that f!62l) is invariant under the following discrete transformations: 

/ - -/ (71) 

t -> -t (72) 

and 

x — > —x (73) 
Taking these reflections into account, we have 

G 4 -> -G A (74) 

G 2 -> -d (75) 

and 

G 3 - -G 3 (76) 
We also note that x — ► — x will make G± — G 2 . 
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d d 
dt dx 


(77) 


'as 


(78) 


^ + ^ + f ^ 

dt dx df 


(79) 


^ +t ^ + Pf ® 
dt dx df 


(80) 


4+< 


(81) 



In order to obtain group invariant solutions of ()62|) (with / = f(t,x) explicitly) we only need 
consider the following symmetry combinations (Msomi 2003): 

Gi : 

G\ + G4 
G3 + PG4 

G! + G 2 + (3G 4 

where (3 is a real arbitrary parameter. All other solutions of (jfi2j) obtained via linear combina- 
tions of the point symmetries ()63J ) - (J66|) can be obtained from the solutions we report here. We 
consider each in turn. 

6.1 Invar iance under G\ 

Using G\ we determine the invariants from the invariant surface condition 

y =f = f (82) 

which yields 

y = x-t (83) 
f = V. (84) 

Now by using the above transformation, in equation ()62|). the partial differential equation is 
reduced to 

k 2 V 2 = (85) 

ic. 

V = 0. (86) 
Thus only trivial travelling wave solutions are possible. 
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6.2 Invar iance under G4 



In this case the system has invariants given by 

y = t (87) 

V = x. (88) 

From the above transformation, we cannot reduce (|62|) . There is no ordinary differential equa- 
tion to solve. 



6.3 Invariance under G\ + G4 

Here, the invariants are 

y = x-t (89) 
/ = VeK (90) 

The partial differential equation is reduced by the transformation to the form 

{k{\ + k)V 2 + V^-V{2KV v + V vu )) = (91) 

which is linearisable and has solution 

where Co and C\ are arbitrary constants, whence 

/ - -p(fi±fl. + ^-fl±fi + g + ^) w 

and so we have "time-boosted" travelling wave solutions. With h given by (joTj) and / by (|93|) 
it is possible to show that there are regions for which \x > and p > 0. This is a desirable 
feature because we expect that barotropic matter in cosmological models should have positive 
pressures and positive energy densities. 

6.4 Invariance under G3 + flG^ 

For this case the system has invariants given by 

y = x 2 -t 2 (94) 
/ = (x + tfV (95) 
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which reduces the partial differential equation to 

k 2 V 2 + A(V y (l + (5)) 2 + SyVyVyy - AVV y k(l + (5)- WV yy ky + WV yy f5(l - (3) = 0. (96) 
This equation is an ordinary differential equation with the single symmetry 

* = v w (97) 



In order to solve the differential equation (J96J) with one symmetry, we first try to reduce the 
order of this equation and see if the resulting equation is easily integrated. 

The differential equation ()96j) has the following reduction variables from (|97j) : 

r = y (98) 

Q = y- (99) 

Using these invariants in (|96j) . we have 

dq _ -k 2 - 4g 2 (l + /3) 2 + 4g(l + /3) _ 2 

dr 8pq-4qkp-4:(-l+p)P) Q ' 1 ' 

In this case, we find that the first order equation cannot be easily integrated. We would 
have to resort to numerical solutions. However, some special solutions can still be found as 
demonstrated later. 



6.5 Invariance under G\ + G2 + fiG^ 
In this final case the invariants are 

y = x (101) 
/ = Ve%. (102) 
The partial differential equation is reduced by this transformation to the form 



2 + 4 A; VVyy - p 2 V 2 - k ({3 2 + 4k) V 2 = (103) 



which has general solution 

<3 2 +4fc / (Ci si 
= 9—5— I -L-t — 

(P 2 + 4k) 



P 2 +4k 

V = 2^r~ -^J — JL — -L. ZjL } (104) 



where C\ and C 2 are arbitrary constants and 



2 kcc , , 



Thus a solution to (Jo^j) is given by 



/3 2 +4fc *£ / ( Cj sin ^ — C 2 cos ■0) ^ s ' 



f = r^**{- ~ • (106) 
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7 A Particular Solution 



In order to investigate the physical properties of the spacetimes obtained, we take a special 
case of the intractable result in ffi.4| that of /3 — 0. In this case the new independent variable 
is given by 

u = t 2 - x 2 (107) 

and the new functional form is 

f(t,x) = f(u) = f(t 2 -x 2 ). (108) 

As a result, (fH2"|) becomes 

d 2 f df 

4m-4 + 2f- + kf = 0. (109) 
air cm 

We now define new variables a and q via 

/(u) = a(«)w 1/4 o = w 1/2 (110) 

and obtain the equation 



dq 2 dq \ \2 
We now distinguish between the two cases k > and k < 0. 

For /c > 0, (jlllj) can be simplified to 



w 2 — + w—+\w 2 -[-\ I o = (112) 



d 2 a da / 2 /l x2 ' 

by the transformation = w ■ Equation (|112|) is the Bessel equation of order | and its solu- 
tions are the linearly independent Bessel functions J±(w) and J_i(w) which may be expressed 
in terms of elementary trigonometric functions in the following way: 



Ji(w) = \ sin w 

2 \ TTW 



2 

J_i(w) = \ COS 10. 

2 \ TIW 



It is a pleasing feature of this model that ()62|) admits Bessel functions as solutions since many 
realistic phenomena are governed by the Bessel equation. The solution to the differential 
equation (|62|). in the original variables t and x, is 



f(x, t) = A sin Jk(t 2 -x 2 ) + B cos Jk(t 2 - x 2 ) (113) 
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where A and B are arbitrary constants. For the solution ()113|) the quantity v is given by 

t 

tanhw = — 

x 

so that the exact solution corresponding to (jll3|) is not conformally flat. 

For k < 0, we make the substitution k = —f3 2 followed by W = (3z. Equation then has 

the form 

, d 2 a da ( 9 /1\ 2 \ 

W 2 + W \W 2 + - ] =0 

dW 2 dW \ \2j J 

which is the modified Bessel equation of order |. This differential equation admits the modified 
Bessel functions 



Ii = J— — sinhW (114) 

2 V TlW 



I_i = J^-coshVy (115) 



as linearly independent solutions. The solution to the field equation ([62)1. in terms of the 
original variables t and x, may be written as 



f(x, t)=A sinh sjk (t 2 -x 2 ) + B cosh yA; (t 2 - x 2 ) (116) 

analogous to ()113|) . 



8 Conclusion 

The notion of conformally mapping a given line element to a new metric, such as the conformally 
related metric (0), is potentially a very fertile avenue in generating new solutions to the Einstein 
field equations. However, the actual number of exact solutions found using this algorithm is 
very low as pointed out by Stephani et al (2003). The Petrov type D spacetime turns out to be 
a rare metric which is amenable to this approach, and we have therefore focussed our attention 
on this line element (first considered by Castejon-Amenedo and Coley (1992)). The Petrov 
type D model investigated is physically well behaved as there is a barotropic equation of state, 
the weak and strong energy conditions are satisfied and it is amendable to a simple two perfect 
fluid interpretation. 

In order to demonstrate a new solution we had to find new functions h and / which define the 
conformal function XJ in (jHJ). Two categories of solutions arise naturally in the analysis. In the 
first category (h = 0) we found a general class of solutions in terms of elementary functions. 
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This contains the solutions of Castejon-Amenedo and Coley (1992) as special cases. In the 
second category (h ^ 0) we invoked the method of Lie group analysis to solve a master field 
equation which allowed us to obtain various exact solutions for the metric under consideration. 
Four cases arise depending on the Lie symmetry generators (we ignore the non-reducible case). 
Three of the cases were readily solvable, in general, in terms of elementary functions. However, 
we were only able to provide special solutions in the remaining case. 

This was analysed further for physical plausibility by considering a particular solution. We 
demonstrated that in this case that the solution can be expressed in terms of Bessel and 
modified Bessel functions. The positivity of pressure and energy density in this case ensured 
its relevance for the description of some cosmological processes. All solutions found via the 
Lie method could be used in both h = and h ^ cases. This treatment demonstrates 
the importance of the method of Lie group analysis in seeking solutions to the Einstein field 
equations. 
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